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Many-body correlations govern a variety of important quantum phenomena such as the emergence 
of superconductivity and magnetism. Understanding quantum many-body systems is thus one of 
the central goals of modern sciences. Here we demonstrate an experimental approach towards this 
goal by utilizing an ultracold Rydberg gas generated with a broadband picosecond laser pulse. We 
follow the ultrafast evolution of its electronic coherence by time-domain Ramsey interferometry 
with attosecond precision. The observed electronic coherence shows an ultrafast oscillation with 
a period of 1 femtosecond, whose phase shift on the attosecond timescale is consistent with many- 
body correlations among Rydberg atoms beyond mean-field approximations. This coherent and 
ultrafast many-body dynamics is actively controlled by tuning the orbital size and population of the 
Rydberg state, as well as the mean atomic distance. Our approach will offer a versatile platform to 
observe and manipulate non-equilibrium dynamics of quantum many-body systems on the ultrafast 
timescale. 


INTRODUCTION 

Atomic, molecular and optical physics with advanced 
laser technologies has recently emerged as a new platform 
to study and possibly simulate quantum many-body sys¬ 
tems m- One of its latest developments is the study 
of long-range interactions among ultracold particles, in 
which dipolar quantum gases HIS], ion crystals [SHE], 
polar molecules [SHII] and Rydberg atoms [HHS] allow 
for revealing the effects of many-body correlations. Ry¬ 
dberg atoms distinguish themselves by their large dipole 
moments and tunability of the strength and nature of 
dipolar interactions [HIS]. 

The active electron in a Rydberg atom moves in a 
macroscopic orbital whose size could range from sub¬ 
micrometre to several tens of micrometres [16j . The re¬ 
sulting large dipole moments yield a strong interaction 
U{r) of either van der Waals or resonant dipolar char¬ 
acter between a pair of Rydberg atoms separated by r. 
This interaction, which can be tuned in various ways 
features the emergence of atomic many-body correlations 
in an ultracold gas. As a prominent example, these in¬ 
teractions shift the energy levels of atoms encompassing 
a given Rydberg atom, so that additional Rydberg ex¬ 
citations are suppressed for U{r) larger than the exci¬ 
tation line width. This effect, referred to as Rydberg 


blockade [niiiin], results in spatial or temporal cor¬ 
relations among the atoms, as recently demonstrated in 
various settings [HI [T8U26] . Its applications include, for 
example, the realization of universal atomic and photonic 
logic gates for quantum information processing [niiiTi- 
m- The Rydberg blockade has also turned out to be 
an outstanding new resource to investigate many-body 
problems [32] • The blockade condition determines the 
smallest distance between neighbouring Rydberg atoms, 
which is typically of the order of several micrometres. 
However, Rydberg excitations can also be induced at 
smaller interatomic distances by, for example, detun¬ 
ing the excitation-laser frequency from an atomic res¬ 
onance [331 - 138] . In this regime, the character of the 
correlations changes, resulting in the facilitated forma¬ 
tion of aggregates consisting of large number of Rydberg 
atoms [5iH55] . 

A complementary approach to the correlations induced 
by the interparticle interactions consists in studying the 
temporal evolution of electronic coherences of the Ryd¬ 
berg atoms. Using a broadband picosecond laser pulse, 
Rydberg excitations can be induced over a wide range 
of interatomic distances from < 1 fim to the isolated 
atom limit. The number of Rydberg excitations per unit 
volume can thus be larger than the ones in the block¬ 
ade regime by two orders of magnitude. Although pi¬ 
cosecond and femtosecond laser pulses have been previ- 
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FIG. 1: Schematic diagram of the experimental setup and Rydberg excitation, (a) Sketch of the experimental setup. 
The dipole-trap laser is turned off 2 /rs before the irradiation of the picosecond pulses, to avoid 2-1-1 multiphoton ionization 
induced by a combination of the picosecond pulses and the trapping laser beam, (b) Two-photon pump (probe) excitation of 
the Rb atom to its Rydberg states, (c) The Rydberg states can be resolved by field ionization with a slowly ramped electric 
field (see Methods section ‘Rydberg excitation and detection’). Here, the Rydberg population and estimated peak atom density 
were 1.2 ± 0.1 % and ~ 4 x 10^° cm“^, respectively (see Methods sections ‘Rydberg excitation and detection’ and ‘Estimation 
of the atom density’ for these population and density estimations). The small peak around 5.6/rs could be assigned to free 
ions generated by population redistribution m and/or direct multiphoton ionization, (d) Sketch of the two-body interaction 
and pulse excitation accompanied by a plot of the nearest-neighbor distribution of Rydberg atoms for the peak atom density 
of n = 1.3 X 10^^ cm“® (pink solid and red dashed traces) and the averaged density over the whole atoms (light-blue solid and 
dark-blue dashed traces) estimated for the present ensemble of the Rb atoms (see Methods section ‘Estimation of the atom 
density’). Here, the pink and light-blue solid traces are obtained by a Monte-Carlo simulation, whereas the red and dark- 
blue dashed traces show analytical results for homogeneous distributions. The difference between the light-blue and dark-blue 
average density traces results from the difference between their Gaussian and homogeneous density distributions, respectively, 
over the whole ensemble. At the peak atom density, the average nearest-neighbour distance is given by 0.5/rm, whereas for the 
averaged density we obtain 0.87 fim. 


ously used to observe electronic wave packets in isolated 
Rydberg atoms [SSI SO] and the dephasing due to two- 
body interactions im, here we exploit them to explore 
how coherent dynamics evolves in a many-body regime. 
The strong interactions in our ultracold Rydberg gas in¬ 
duce an electronic dephasing on the picosecond timescale, 
which is directly observed in a time-domain Ramsey in- 
terferogram oscillating with a period of ^ 1 fs and phase- 
shifted on the attosecond timescale |30| by the Rydberg 
interactions. We measure this minute phase shift and 
the dephasing directly, and compare them with theoreti¬ 
cal simulations based on nearest-neighbour interactions, 
a mean-field model and many-body correlations to reveal 
effects indicating atomic correlations beyond a mean-field 
description. The two-body interaction energy at around 
1 ^m distance, which can be easily accessed in our ap¬ 
proach, exceeds the average kinetic energy by many or¬ 
ders of magnitude. This regime compares favorably with 


previous experiments in thermal vapour cell experiments 
with nanosecond laser pulses [351 US] ? where the nearest- 
neighbour distance is also < 1 /rm. In thermal cells, the 
Rydberg interaction energy is comparable to the aver¬ 
age kinetic energy, so that the measurement of coherent 
evolution of the interaction dynamics is strongly affected 
by the thermal atomic motion, as will be discussed later 
quantitatively. 


RESULTS 
Experimental set-up 

Figure & shows the schematics of our experimen¬ 
tal setup. A cold ensemble of ®^Rb atoms is pre¬ 
pared in an optical dipole trap with temperature and 
highest peak atom-density estimated to be ^70 p,K and 
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~1.3xl0^^ cm“^, respectively (see Methods sections ‘Es¬ 
timation of the temperature’ and ‘Estimation of the atom 
density’ for these temperature and density estimations). 
On the timescales of interest, atomic motion can be ig¬ 
nored in the frozen gas regime [HI 03] . The atoms are 
optically pumped to the hyperfine state F=2, mF=-l-2 
of the ground state 5Si/2 and excited to Rydberg states 
via a two-photon transition using broadband picosecond 
laser pulses with their centre wavelengths tuned to ~ 779 
and ^ 481 nm (Fig. [^), hereafter referred to as the in¬ 
frared and blue pulses, respectively. The dipole-trap laser 
is turned off 2 /is before the irradiation of the infrared and 
blue pulses to avoid 2-1-1 multiphoton ionization induced 
by a combination of the infrared and blue pulses with 
the trapping laser light. The infrared and blue pulses, 
and the optical pumping beam are circularly polarized 
in the same direction with respect to the magnetic held, 
suppressing excitations to the S Rydberg states, so that 
the state 1 ^ 05 / 2 , mj=-1-5/2 is mostly populated, where v 
is a principal quantum number. The population of the 
Rydberg state is measured by held ionization [TB]. De¬ 
tails of the held ionization are described in Methods sec¬ 
tion ‘Rydberg excitation and detection’. The maximum 
population of the i^D states in the present experiment is 
not > 5 %, to suppress photoionization by the picosecond 
laser pulses (see Supplementary Note 1 for the effects of 
ions). More details on the atom preparation, the esti¬ 
mation of the atom density and the temperature, as well 
as the Rydberg excitation are described in Methods and 
Supplementary Note 2. 

The bandwidth of our Rydberg excitation is deter¬ 
mined from a held-ionization spectrum exemplihed in 
Fig. (see Supplementary Note 3), in which the excita¬ 
tion is tuned to the 42 D 5/2 state, which is the main tar¬ 
get state of the current experiment. The held-ionization 
voltage is ramped up slowly enough on the 5 /iS timescale 
to resolve neighbouring Rydberg levels in such held- 
ionization spectra. The bandwidth of the excitation is 
~ 150 GHz (full width at half maximum, FWHM) and is 
much larger than those of the continuous-wave, nanosec¬ 
ond, and sub-nanosecond pulsed lasers employed in pre¬ 
vious ultracold Rydberg experiments [m 03114111461147]. 
The held-ionization spectrum indicates that the band¬ 
width of the Rydberg excitation is larger than the energy 
separation of neighbouring Rydberg states and is accord¬ 
ingly wide enough to remove the Rydberg blockade. As 
schematically shown in the upper panel of Fig. [^, our 
picosecond laser pulses can excite a pair of Rb atoms 
simultaneously to the Rydberg states ly = 42 even at 
interatomic distances shorter than 1 /im. 

Observation of many-body electron dynamics 

The interaction among the Rydberg atoms is observed 
by time-domain Ramsey interferometry with a pair of the 


Rydberg excitations, hereafter referred to as the ‘pump’ 
and ‘probe’ excitations, whose delay was stabilized on the 
attosecond timescale with our homemade optical inter¬ 
ferometer 05] (see Methods section ‘Time-domain Ram¬ 
sey interferometry’). By scanning the delay time r be¬ 
tween the pump and probe excitations on the attosecond 
timescale, we measure the Ramsey oscillation of the pop¬ 
ulation integrated over all Rydberg states, which remains 
after the probe excitation, by the field ionization. The 
field-ionization voltage in these Ramsey measurements 
is ramped up rapidly enough on the 100 ns timescale to 
avoid subsequent decay processes of the Rydberg states. 
For the excitation tuned to the 42 D 5/2 state as shown in 
Fig. a single Rydberg state \ v) with 1 / = 42 is predom¬ 
inantly populated. The Rydberg population thus mea¬ 
sured oscillates as a function of r with a frequency close 
to the transition frequency between the 5S and Rydberg 
states 132 . This is in contrast to the standard Ramsey 
interferometry, in which the signal oscillates with a fre¬ 
quency close to the detuning frequency of the excitation 
laser from the atomic transition \ni I23j . That is, in the 
absence of interactions, the population in the Rydberg 
state Pv{t) is given by 

Pu{t) oc 1 -I- cos{E^T/h) (1) 

and oscillates with the frequency E^/h^ where E^, is 
the energy of the Rydberg state \v) measured from the 
ground state 5S (see refs 42, 49 and Supplementary Note 
2), and h is the Planck constant h divided by 27r. This 
oscillation is identical to the temporal oscillation of the 
Rydberg state |^), except that the real-time t is replaced 
by the pump-probe delay r. Therefore, the Ramsey oscil¬ 
lation for 1 / = 42 corresponds essentially to the temporal 
oscillation of the Rydberg wave function |z/ = 42) and 
to the recurrence motion of an electronic wave packet, 
which is composed of the 5S and Rydberg state ^ = 42 
superposed coherently by the excitation pulses. Here we 
investigate how the Rydberg interactions affect this co¬ 
herent electron dynamics. 

In a simplihed mean-field approach, the Rydberg in¬ 
teractions change only the energy of the Rydberg state 
|z/) and accordingly the oscillation period of Pv{t). This 
results in a phase-shift of the Ramsey oscillation accumu¬ 
lated in the delay time r. In addition, as the atoms are 
randomly distributed in the ensemble, their energy levels 
El, are shifted randomly by the interactions, making the 
periods of their Ramsey oscillations different from each 
other. Therefore, the measured signal is the superposi¬ 
tion of many oscillations with different periods, so that 
its contrast is expected to decay as a function of r due 
to Rydberg interactions. 

Figure [^-c show examples of the Ramsey oscillations 
for the 42 D 5/2 state with a population of 3.3 ± 0.1 %. As 
mentioned above, these oscillations correspond to the ul¬ 
trafast recurrence motion of the electronic wave packet 
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FIG. 2: Ramsey oscillation for the 42D state, (a - c) The field ionization signals for the higher-density (red traces) and 
the lower-density (blue traces) ensembles are plotted as functions of the pump-probe delay r scanned over a range of ~ 3fs 
around 75, 245 and 500 ps, showing clear oscillations. The signal intensities are normalized by the mean value of the sinusoidal 
function fitted to each oscillation. The origin of the pump-probe delay r = 0 is arbitrary and is taken to be the left edge of each 
figure, (d) The contrasts of these oscillations are plotted as functions of r. The oscillatory structures as functions of r could be 
partly attributable to the recurrence motion of a wave packet composed of the 42D and its neighbouring Rydberg states (see the 
field-ionization spectrum shown in Fig. Ic). The shortest period of this recurrence motion is evaluated from their level spacing 
to be ~ 10 ps, which is not resolved in these plots. More details on the oscillatory structures are described in Supplementary 
Note 4. The contrast decays clearly as a function of r for the higher-density ensemble, whereas for the lower-density one a 
decay is not visible, (e) The phase shift of the higher-density ensemble from the lower-density one is plotted as a function of r. 
The oscillatory structure as a function of r is partly attributable to the wave-packet motion. The error bars represent the s.d. 


with a period of ^ 1 fs. For each oscillation, we mea¬ 
sured the field-ionization signals of two atomic ensembles 
with different peak densities estimated to be ^ 1.3 x 10^^ 
and ^ 4 X 10^° cm“^ alternately to suppress systematic 
uncertainties, scanning r in steps of ~ 30 as (see Meth¬ 
ods section ‘Estimation of the atom density’ for these 
density estimations). We obtained the contrasts and 
phases of the measured oscillations by sinusoidal fitting, 
as shown in Fig. [^-c (see Methods section ‘Time-domain 
Ramsey interferometry’). Figure]^ shows that the con¬ 
trast is approximately constant for r up to ~ 500 ps 
for the lower-density ensemble (blue-circle data points), 
indicating that the effects of the interactions are neg¬ 
ligibly small. This result is consistent with the inter¬ 
action strength estimated from the present atom den¬ 


sity and the two-atom potential curve presented in Sup¬ 
plementary Fig. 1. Hereafter, we take these contrasts 
and phases measured in the lower-density ensemble as 
references to be compared with those measured in the 
higher-density ensemble. The oscillatory structures on 
the ^ 10 ps timescale at r ~ 130 — 170 ps seen in Fig. 2d 
could be partly due to the recurrence motion of a wave- 
packet composed of the Rydberg state u = 42 and the 
traces of its neighbouring states seen in Fig. Ic (see Sup¬ 
plementary Note 4 for details). 

Figure [2ji shows that the contrast decays as a func¬ 
tion of T for the higher-density ensemble (red-circle data 
points). The phase shift of the higher-density ensemble 
from the lower-density one also changes as a function of 
T as seen in Fig. The offset of this phase shift at 
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T = 0 is essentially due to the difference between AC- 
Stark shifts of the atomic levels in the higher- and lower- 
density ensembles. Slight differences between the sizes, 
shapes and positions of the higher- and lower-density en¬ 
sembles could lead to their different AC-Stark shifts (see 
Supplementary Note 2 for more details on the origin of 
this zero-delay offset of the phase shift). 

Figure shows the results of measurements of the 
contrast ratio between the two different densities given 
above, hereafter referred to as ‘Ramsey contrast’, and 
the phase shift for the 42 D 5/2 state. These results are 
plotted as functions of r for the two different Rydberg 
populations Pe of 1.2±0.1% and 3.3±0.1%, showing that 
the Ramsey contrast decays, and the phase shift is accu¬ 
mulated as a function of r. 

It is evident from Fig. that the contrast decay and 
the phase shift are enhanced when the Rydberg popula¬ 
tion is increased from ^1.2 % to ~3.3 %. The strength of 
the interactions is tuned by varying the Rydberg quan¬ 
tum number ly and the atom density. Figure shows 
Ramsey contrasts as functions of r for three different Ry¬ 
dberg levels iz = 38, 42 and 50. It is seen from this figure 
that the contrast decay is accelerated by increasing the 
principal quantum number v of the Rydberg level. The 
dependence of the Ramsey contrast decay on the atom 
density is shown in Fig.[^, in which the contrast decay is 
accelerated by increasing the atom density (see Supple¬ 
mentary Note 5 for the estimation of the atom densities 
plotted in the abscissa of this figure). From these com¬ 
bined measurements as functions of the Rydberg popu¬ 
lation, principal quantum number and atom density, we 
conclude that the observed contrast decay and the phase 
shift are induced by Rydberg interactions. 

The origin of the observed behaviour of the contrast 
decay and the phase shift is further investigated in Fig.[^ 
Here we compare the experimental data with the Ram¬ 
sey contrast decays and the phase shifts calculated for 
nearest-neighbour interactions without considering inter¬ 
actions among three or more Rydberg atoms (solid lines). 
The zero-delay offset of the calculated phase shift is ar¬ 
bitrary and adjusted so that the average of the first ten 
data points is equal to the calculated phase shift averaged 
over the delay window for those ten data points. Results 
for nearest-neighbour interactions are obtained by simple 
numerical calculations following those performed in pre¬ 
vious Ramsey studies [23111], in which the time-domain 
Ramsey oscillations modulated by nearest-neighbour in¬ 
teractions are averaged over the distribution of nearest- 
neighbour distances (see Supplementary Note 6 for de¬ 
tails). Here, the Ramsey oscillation is obtained by solv¬ 
ing a Schrodinger equation with a Hamiltonian for two 
interacting atoms. We assume a van der Waals interac¬ 
tion of the form U{r) = —C^jr^ with the coefficient Cq 
being an adjusting parameter. It is seen in the insets of 
Fig- Hi ,b that the calculated Ramsey contrast decay be¬ 
comes slightly faster as the Cq value is increased, but does 


not grow beyond the decay thresholds 0.988 and 0.967 
determined by the Rydberg population pe ~ 1.2% and 
3.3%, respectively (for details on these thresholds, see 
Supplementary Note 6). Similary, the calculated phase 
shift is almost constant as the Cq value is increased, as 
shown in Fig. & ,d. These features are not altered by 
changing the character of the interaction such as van der 
Waals and dipole-dipole (see Supplementary Note 6 for 
details). The contrast decays and phase shifts observed 
experimentally are thus clearly larger than expected for 
nearest-neighbour interactions. It is therefore concluded 
that our experimental observation of the Ramsey con¬ 
trast decay and phase shift on the attosecond timescale 
demonstrates interactions among more than two Rydberg 
atoms, and the effect of those interactions on the electron 
dynamics can be actively controlled by tuning the pop¬ 
ulation and principal quantum number of the Rydberg 
level, as well as the atom density. 

Test of a mean-field model 

To model these observations, we first apply a mean- 
field model (see Supplementary Note 7 for its details). 
As in the calculations with nearest-neighbour interac¬ 
tions above, we assume a van der Waals interaction of 
the form U{r) = —C^jr^ with the coefficient Cg being 
the only fitting parameter, which is optimized to repro¬ 
duce the measured Ramsey contrast in Fig. Wp by a least- 
squares fitting. The outline of the least-squares fitting is 
given in Supplementary Note 8. The results of the mean- 
field simulations are plotted in Fig. [^(dashed lines). Fig¬ 
ure ,b show that reasonable agreements can be found 
for the Ramsey contrast between the measurements and 
the mean-field simulations with Cg = 1.9 GHz/rm®. This 
value for Cg yields the phase shifts simulated by the 
mean-field model in Fig. ,d, where the zero-delay offset 
is adjusted in the same way as in the calculations with 
nearest-neighbour interactions above. Figure shows 
that the mean-field simulation yields a phase-shift larger 
than the measured one by a factor of ~ 4, failing to 
reproduce our observations. This discrepancy between 
measured and simulated phase shifts is not improved by 
assuming a dipole-dipole interaction, a hybrid form of 
a dipole-dipole and a van der Waals interaction or by 
introducing anisotropic interactions (see Supplementary 
Figs 2 and 3). The closer agreement between the mea¬ 
sured and simulated phase shifts for pe ~ 3.3 % in Fig. 3d 
than for pg ^ 1.2% in Fig. 3c is understood as follows. 
The Gaussian atom-density distribution of our experi¬ 
mental setup leads to the saturation of the phase shift 
at longer pump-probe delays. This is because the rapid 
decrease of the atom density in the Gaussian tails results 
in the suppression of the contribution to the phase shift 
from atoms distant from the centre, and therefore the 
phase shift does not grow afterwards. This saturation 
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FIG. 3: Measured Ramsey contrasts and phase shifts compared with those expected for nearest-neighbour 
interactions and mean-field approximations. The black diamond-shaped and blue circle data points show the measured 
Ramsey contrasts (a and b) and the phase shifts (c and d) for the popnlation pe ~ 1.2% (a and c) and ~ 3.3% (b and 
d) of the 42 D 5/2 state, respectively. In a and b, the measured Ramsey contrasts are compared with the simulated ones for 
nearest-neighbour interactions with a Ca coefficient of 2 GHz pm® (green and yellow solid lines) and 50 GHz pm® (purple and 
orange solid lines), respectively Results of the mean-field simulations are presented by magenta and red dashed lines in a 
and b, respectively Similarly, the measured and simulated phase shifts are compared in c and d. The zero-delay offset of the 
simulated phase shift is arbitrary and adjusted so that the average of the first ten data points is equal to the calculated phase 
shift averaged over the delay window for those ten data points. In the insets in a and b, the Ramsey contrasts simulated for 
the nearest-neighbour interactions (solid lines) are vertically magnified to show their convergence as the Ge value is increased. 
Grey dot-and-dash lines show their population-dependent thresholds given by 1 — Pe- In the insets in c and d, the phase 
shifts simulated for the nearest-neighbour interactions (solid lines) are magnified vertically The peak atom density is set to 
~ 1.3 X 10^^ cm“® in these simulations. The error bars represent the s.d. 


is reached both by the measured and simulated phase 
shifts for Pe ~ 3.3 % within our measurement time 500 ps 
(unlike the case with pe ^ 1.2 % and therefore weaker in¬ 
teractions), giving their closer agreement for pe 3.3% 
than for Pe ~ 1.2%. This saturation effect is explained 
more quantitatively in Supplementary Note 9. 


Beyond mean-field analysis 

Next, we apply an exactly solvable theory model [501- 
151] to the observations. Details of this model are pre¬ 
sented in ref. 54. Briefly, we represent each atom as a 
two-level system, which is a pseudo-spin system, consist¬ 
ing of a ground state |g) (= |4,)) and an excited Rydberg 
state |e) (= |t)) with energies iHg and E^., respectively. 
The experiment consists of four stages: (i) the pump ex¬ 
citation, (ii) an evolution with an A^-atom Hamiltonian 
H, (iii) the probe excitation and (iv) the population mea- 
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FIG. 4: The principal quantum number and atom density dependences of the Ramsey contrast, (a) Measured 
Ramsey contrasts are plotted as functions of r for three different Rydberg levels v = 38,42 and 50. The estimated populations 
Pe and peak atom densities are Pe ~ 3.2 % and ~ 1.2 x 10^^ cm“^ for v = 38, Pe ~ 3.3 % and ~ 1.3 x 10^^ cm“^ for v = 42, and 
Pe ~ 3.1 % and ~ 1.2 x 10^^ cm“^ for u = 50, respectively (see Methods sections ‘Estimation of the atom density’ and ‘Rydberg 
excitation and detection’ for these population and density estimations). The simulations indicated by the black, blue and red 
solid lines yield adjusting parameters to be Ce = 8, 34 and 103 GHz pm® for v — 38,42 and 50, respectively. The interaction 
strength in these simulations is limited below 75 GHz, which is the half width half maximum of the pump excitation, and the 
peak atom density is set to the estimated density for each Rydberg level in these simulations. It should be noted that several 
Rydberg states are excited for iz = 50 (see Supplementary Fig. 4). However, we consider an excitation to a single Rydberg 
state, to perform the simulations for all of the three Rydberg levels. Each error bar in the inset represents the average over 
the error bars (the s.d.) of all data points for each Rydberg level, (b) Ramsey contrast is measured as a function of the peak 
atom density at two different pump-probe delays r = 300 and 510ps for v = 42 with its population being ~ 3.5 %. The vertical 
error bars represent the s.d. and the horizontal error bars arise from the density calibration (see Supplementary Note 5 for the 
estimation of the atom densities plotted in the abscissa). 


surement by field ionization. The A'^-atom Hamiltonian 
in stage (ii) includes the atomic energies and the inter¬ 
actions U{rjk) between atoms in the Rydberg states as 
follows: 


^ ] 
H = ^ huj - 

j=i 


N-l N 

k—1 j>k 


1 + at 


,( 2 ) 


with oj = {Eg — Eg)/h the atomic-resonance frequency 
and (t| the Pauli matrix, which is used to represent the 
internal states of the pseudo-spin at position j (see Meth¬ 
ods section ‘Outline of the exactly solvable model simu¬ 
lation’ for more details) [50H54] . 

An exact solution for the time evolution with a Hamil¬ 
tonian of the form of equation ([^ has recently been 
presented in refs 50-54 (see also similar numerical ap¬ 
proaches in refs 55,56). This allows for deriving an ex¬ 
pression for the exact time evolution of the Ramsey sig¬ 
nal P{t) for any strength of interactions [5DH51] . For 
any given atom j interacting with — 1 neighbouring 


Rydberg atoms, one obtaines 


Pjir) 


2pgPe^< 1 + e 


2(cjr+0) 


N 


k^l 


(3) 


where Pg and Pe are the ground- and Rydberg-state pop¬ 
ulations, respectively, produced by the initial pump ex¬ 
citation, Ajk = U{rjk)/fi describes a frequency shift in¬ 
duced by the interaction between atoms j and fc, and (j) 
is the phase offset arising from the AC-Stark shifts dur¬ 
ing the picosecond pulse excitations (see Supplementary 
Note 2). It is seen from this equation that this model con¬ 
siders different clusters of interactions superposed coher¬ 
ently instead of averaged as performed in our preceding 
mean-held analysis. This coherent superposition leads to 
correlations among different atoms, as exemplihed for a 
two-atom correlation in Fig. (ref. 54). 

Further analytical progress is possible by using an ap¬ 
proximation, hereafter referred to as a ‘continuum ap¬ 
proximation’, in which a continuum function n(r) is con¬ 
sidered for the density distribution of Rydberg atoms (see 
Methods section ‘Continuum approximation’ for details). 
Briehy, we assume the density distribution to be homo- 
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FIG. 5: Two-atom correlation simulated by the exact 

model. The two-particle correlation function is calculated for 
two atoms in the Gaussian atom density distribution of the 
present Ramsey measurements for the Rydberg level i/ = 42. 
The one atom is located at the centre of the distribution, and 
the distance to the other atom is set to 1 pm (a) and 1.3 pm 
(b), respectively. The Ce coefficient, the Rydberg population 
and the peak atom density are set to 34 GHz pm®, ~ 3.3 % and 
~ 1.3 X 10^^cm“®, respectively. The interaction strength is 
limited below 75 GHz, which is the half width half maximum 
of the pump excitation. The coherent superposition of differ¬ 
ent clusters of interactions leads to dephasing and therefore 
the global decay. 


geneous in a small volume around a particular position r 
(ref. 54). This approximation leads to the following ex¬ 
pression for the Ramsey signal averaged over the whole 
atomic ensemble 


P{t) = 2pgpe {1 + | 5 (r)| cos(a;r -b a(r) -b (p)} . (4) 


The contrast decay |p(t)| and phase shift a(r) are then 
obtained from p(r) = |p(T)|e*“*^’’\ which, in the case of an 
isotropic van der Waals interaction and for r > 27r/wB, 
is given by 


9{t) = 


/TTrin 


dn\lln ) 


„-crTi(l-dVr 



(5) 

Here, wb = 27r x 75 GHz is the half width half maxi¬ 
mum of the pump excitation, Up is the peak atom den¬ 
sity, the coefficient 2-\/ln (up/n)/( i/ttUp) results from the 
Gaussian atomic density distribution in the experiment 
and (T„ = Pen^/SnCe/h (tt/S) is the decay constant [M] . 
Equation ([^ predicts that the Ramsey contrast decays 


approximately as a stretched exponential e~°‘^ with 
a square-root dependence on the pump-probe delay r. 
This square-root dependence is characteristic of the van 
der Waals interaction [^. The decay constant reads 
a where riav = 2~^/^np is the average density for 

a Gaussian distribution. 

Figure shows comparisons between (i) the experi¬ 
mental data for the contrast decay and phase shift as 
functions of r for the 42 D 5/2 state and (ii) their nu¬ 
merical results based on the analytical continuum ap¬ 
proximation equation Q (solid lines). The curves for 
the exact numerical results agree well with the analytical 
ones for 00 , so that they are indistinguishable on 

the present scale of the figure. The figure demonstrates 
that both the exact and analytical results obtained for 
Cq = 34 GHz ^m® agree well with the measured Ramsey 
contrasts and phase shifts for both of the Rydberg pop¬ 
ulations Pe ~1.2% and 3.3% used in the measurement. 
Figure also shows numerical results with the contin¬ 
uum approximation for the Ramsey contrasts and phase 
shifts obtained for a hnite average number of interacting 
Rydberg atoms N = 20 and 40 (semi-transparent lines; 
see Methods section ‘Gontinuum approximation’). The 
agreement between the numerical and experimental re¬ 
sults improves monotonically with increasing N. Figurej^ 
shows the Ramsey contrast and phase shift measured (the 
dark-grey solid lines accompanied by light-grey shaded 
areas) and simulated with the continuum approximation 
(the blue solid lines) at the pump-probe delay r ~ 500 ps. 
It is seen from Fig. that the blue solid line crosses the 
upper boundary of the light-grey shaded area, which is 
a confidence interval of the measured value, at an av¬ 
erage number of interacting atoms around 40, showing 
that more than ^ 40 atoms correlated are necessary to 
reproduce the measured phase shift. 

In the theory model above, the Cq coefficient of the 
van der Waals interaction serves as the only fitting pa¬ 
rameter for the Ramsey contrast in Fig. Wp- The Cq value 
thus obtained is used for calculating the contrast decay 
in Fig. and the phase shifts in Fig. |^,d. Then, the 
zero-delay offset of the phase shift serves as the only ad¬ 
justable parameter for the phase shift in each of Fig. |^,d. 
For distances shorter than the average interatomic sep¬ 
aration, the bandwidth of our Rydberg excitation with 
the picosecond infrared and blue pulses covers multiple 
adiabatic interaction potentials that can hybridize with 
the one correlating asymptotically to the 42D-I-42D limit 
(see Supplementary Fig. 1). Although in principle this 
could affect the contrast decay, the present model based 
on a single effective potential captures the observed dy¬ 
namics both for the Ramsey contrast and the phase shift 
(see Supplementary Note 10 for the effective treatment 
of the interaction potentials). 

We note that the predicted decay constant for the 
Ramsey contrast for pe ^ 3.3 % is approximately given by 
~ X 1.45 X 10^ (Hz)^'^^. The pump-probe delay 
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FIG. 6: Measured Ramsey contrasts and phase shifts compared with the numerical ones based on the analytical 
continuum approximation equation The black diamond-shaped and blue circle data points show the measured 

Ramsey contrasts (a and b) and the phase shifts (c and d) for the population of the A2T>^/2 state being ~ 1.2% (a and c) 
and ~ 3.3% (b and d), respectively. In a and b, the measured Ramsey contrasts are compared with the numerical ones based 
on the analytical continuum approximation for the average number of Rydberg atoms = 20 and = 40 (purple and blue 
semi-transparent lines, respectively) as well as for the limit N ^ oo (black and blue solid lines). The curves for the exact 
solution agree well with the ones for A'^ —>■ oo in the continuum approximation, so that they are indistinguishable on the present 
scale of the figure. The orange and pink shaded areas correspond to 2 s.d. of the Cg coefficient obtained in the least-squares 
fitting. In the insets in a and b, the analytical results in the limit are magnified to show their early quadratic decays at small 
pump-probe delays originating from the limited bandwidth given by 75 GHz, which is the half width half maximum of the 
pump excitation. Similarly, the measured and numerical phase-shifts are compared in c and d. The peak atom density is set 
to ~ 1.3 X 10^^ cm“^ in these simulations. The error bars represent the s.d. 


to reach e“^(~ 0.37) of the initial contrast is thus given 
by ~ 760 ps, which agrees well with our measure¬ 

ment time of 500 ps giving the contrast reduction to 45 %. 

We have further employed ^(t) in equation (i to re¬ 
produce the measured Ramsey contrasts shown in Fig.|^ 
for three Rydberg levels v = 38, 42 and 50, finding good 
agreement. The phase shifts for these levels have also 
been measured and analyzed in Supplementary Note 5. 
It should be noted that several Rydberg states are ex¬ 
cited for 1 / = 50 (see Supplementary Fig. 4). However, 
the simulation with a single Rydberg state gives good 
agreement with the experimental results as shown by the 


solid lines in Fig. 

We have also performed additional calculations of 
the measured Ramsey contrasts and phase shifts us¬ 
ing dipole-dipole interactions, a hybrid form of a 
dipole-dipole and a van der Waals interaction, and an 
anisotropic van der Waals interaction. Although all re¬ 
sults are in qualitative agreement with the experimental 
data, we find that the results for pure van der Waals 
interactions presented above reproduce the experimen¬ 
tal data well with only a single fitting parameter. The 
results of these additional calculations are given in Sup¬ 
plementary Note 10. 
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FIG. 7: Convergence of the simulated Ramsey con¬ 
trast and phase shift as functions of an average num¬ 
ber of interacting atoms. The Ramsey contrast (a) and 
phase shift (b) at r = 500 ps are simulated by the theory 
model with the continuum approximation and are plotted as 
functions of the cutoff radius ro (the lower abscissa) and of 
an average number of interacting atoms within the volume 
V = ^(ro — rg) (the higher abscissa), where rg is the block¬ 
ade radius (see Methods section ‘Continuum approximation’ 
for more details of ro and re). The population of the 42 D 5/2 is 
set to ~ 3.3% in these simulations. The interaction strength 
is limited below 75 GHz, which is the half width half maxi¬ 
mum of the pump excitation, and the peak atom density is 
set to ~1.3x 10^^ cm“® in these simulations. The results with 
the van der Waals interaction are displayed by the blue solid 
lines. The dark-grey solid lines represent the measured Ram¬ 
sey contrast and the phase shift, each of which is the average 
over eight points around r = 500 ps in Fig. [^,d. The light- 
grey shaded area represents 1 s.d. of the average over those 
eight measured values. Similar results with a dipole-dipole 
interaction and a hybrid form of a dipole-dipole and a van 
der Waals interaction are shown in Supplementary Fig. 5. 


Origin of the failure of a mean-field model 


The overestimation of the phase shift by the mean-field 
approximation in Fig. n is intuitively understood as fol¬ 
lows. In our actual measurement, the Ramsey contrast 
of each atom decays due to many different frequencies 
corresponding to different clusters of interactions super¬ 
posed coherently, as is expected from equation (3) of the 


exact theory model. In the mean-field model, however, 
the Ramsey oscillation of each atom has its own mean 
frequency, interacting with ‘the other atoms’ as a whole, 
so that its contrast does not decay, but the decay arises 
only from an ensemble average of Ramsey oscillations of 
many atoms phase shifted from each other. Accordingly, 
those phase shifts need to be overestimated, as seen in 
Fig.|^, to reproduce the measured contrast decay. 

Figure shows a comparison between the phase shifts 
calculated by the mean field and exact models as func¬ 
tions of the Rydberg population at three different 
pump-probe delays t = 20, 50 and 70 ps. It is seen from 
this figure that the mean-field phase shift is larger than 
the exact one for pe < 0.5. It is also seen from this fig¬ 
ure that this difference becomes larger as the effective 
number of interacting atoms becomes larger at longer 
pump-probe delays. The mean-field model thus overesti¬ 
mates the phase shift for the present pe ~ 0.01 and 0.03 
(< 0.5), and r ~ 50 - 500ps, as is intuitively understood 
as described in the preceding paragraph. 


DISCUSSION 

The experiments presented here bear similarities with 
recent investigations on coherent spin-exchange dynam¬ 
ics with rotational states of polar molecules trapped in 
an optical lattice liiin]. In these experiments, the effect 
of single-particle decoherence, which proceeds faster than 
the interaction timescale (~ 10 ms), has successfully been 
circumvented by applying spin echo techniques. In our 
experiment, on the other hand, the interaction timescale 
~ 1 ns is about three orders of magnitude shorter than 
the timescale of single-particle decoherence, which is in¬ 
duced mainly by Doppler broadening and is ^ 1 p,s at 
the temperature ^100 pK estimated for our Rydberg gas 
(see Methods section ‘Estimation of the temperature’). 
The combination of ultrafast and ultracold approaches 
thus provides an effective pathway for isolating the ob¬ 
servation and control of coherent dynamics of a system 
from its single-particle decoherence processes. 

Experiments similar to the ones demonstrated here 
could in principle be implemented in thermal cells |38ll43j 
if one can compensate for the relevant Doppler shift, 
which causes a phase shift comparable to the interaction- 
induced phase shift, and can also compensate the effects 
of atomic motions. In the thermal cell at a room tem¬ 
perature, the atoms are not frozen, but move with ki¬ 
netic energies at that temperature, so that the motion 
would lead to additional phase shifts on the order of 
2ttvtIX ^ 27r X 0.3 under the conditions (A ^ 297 nm, 
T ^ 500ps, V ~ 170 m/s), where A is the wavelength that 
corresponds to the energy difference between the ground 
and Rydberg states. In our current ultracold measure¬ 
ments, on the other hand, the atoms move only 40 pm 
on average during the same duration, and this distance 
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FIG. 8: Comparison between the phase shifts obtained by the mean-field and exact calculations. In a - c, the 

phase shifts obtained by the mean-field (red solid lines) and exact (bine solid lines) calculations are plotted as functions of 
the Rydberg population pe at the pump-probe delays r = 20, 50 and 70 ps, respectively. The frequency shifts induced by 
the interactions are set to [/ = 6.96 GHz, common to all the pairs of atoms. This value corresponds to the energy shift with 
Ce = 34 GHz /rm® at the average internuclear distance in a Gaussian ensemble with its peak density 1.3 x 10^^ cm“®. 


is shorter than A by four orders of magnitude. This dis¬ 
tance is also shorter than the average nearest-neighbour 
distance in our Rydberg gas by four orders of magnitude 
(see Methods section ‘Estimation of the atom density’). 
Our Rydberg gas is therefore safely regarded as a frozen 
gas, so that the effects of the atomic motions and col¬ 
lisions on the measurements of coherent dynamics are 
negligible. 

In our present study, a combination of the contrast and 
phase measurements serves as a useful tool to observe 
the effects of many-body correlations. The correlations 
would be further verified by additional measurements of 
variance in Ramsey signals, as is performed in refs 8,57, 
to observe similar many-body correlation effects. Com¬ 
bining a microscope [niiiii I21j with our experimental 
setup offers another future possibility to observe many- 
body correlations more directly in a spatially resolved 
manner. 

We anticipate promising future applications of ap¬ 
plying ultrashort coherent laser pulses to ultracold Ry¬ 
dberg gases. Using alternative excitation schemes, 
one may investigate beyond mean-field effects in Ram¬ 
sey experiments with more complex Hamiltonians such 
as Heisenberg-type Hamiltonians of interest for polar 
molecules [9l [10] and atomic clocks m in optical lat¬ 
tices. Another application could be the investigation of 
a scenario in which Rydberg electronic wave functions 
are spatially overlapped between neighbouring Rydberg 
atoms [58] . This could lead to new exotic phases in which 
the Rydberg electrons are shared among many nuclei, 
and exchange interactions play key roles in their dynam¬ 
ical properties on the ultrafast timescales. Such a metal¬ 
like many-body Rydberg state would naturally lead to 


Penning ionization quite rapidly. However, Jaksch and 
colleagues [59] have theoretically estimated the lifetime 
of such a metal-like Rydberg state of two ®®Rb atoms 
(iz = 50) to be ~ 100 ns. This is longer than our measure¬ 
ment timescale by more than two orders of magnitude, 
rendering the observation of this state possible using our 
time-domain approach. 


METHODS 
Atom preparation 

A magneto-optical trap (MOT) of ®^Rb atoms was 
loaded from background vapour for 1.4 s. During the sub¬ 
sequent MOT compression for 30 ms, an optical dipole 
trap was turned on. The dipole trap was composed 
of a single l,064nm beam with its power and beam 
waist being ^ 4W and ~ 30^m (1/e^radius), respec¬ 
tively. Polarization gradient cooling was performed for 
100 ms. The trapped atoms were then transferred into 
the F = 1 ground state by switching off the MOT re¬ 
pump laser. After that, we turned off the MOT trapping 
beams and the magnetic field. While keeping the inten¬ 
sity of the dipole trap laser, plain evaporative cooling 
was carried out for 50 ms. During the evaporation pro¬ 
cess, a 76 fiT homogeneous magnetic field was turned on, 
pointing along the direction of the dipole trap laser. For 
the next 200 /rs, the atoms were optically pumped to the 
|5 Si/2 , F = 2, mp = +2) state by using the MOT repump 
beam and a cr'*' beam, which is resonant to the transi¬ 
tion from |5 Si/2 ,F = 2) to |5P3/2,F' = 2) and counter- 
propagates with the dipole trap laser. The dipole-trap 
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laser was turned off 2 /is before the irradiation of the 
picosecond infrared and blue pulses to avoid 2+1 mul¬ 
tiphoton ionization induced by a combination of the pi¬ 
cosecond pulses and the trapping laser beam, whereas the 
homogeneous magnetic field remained on. The picosec¬ 
ond infrared and blue pulses at ~ 779 and ~ 481 nm, 
which propagated collinearly with the dipole trap beam, 
had cross-sections with FWHMs of ~ 130 (100) ^m and 
30 (30) /rm along the x (y) direction (Fig. la). 

Estimation of the atom density 

At first, the atom density in the Ramsey measurements 
was estimated solely from the total number of atoms and 
the size of the atomic ensemble obtained by in-situ ab¬ 
sorption imaging with a CCD (charge-coupled device) 
camera without expanding the atomic ensemble. How¬ 
ever, in contrast to the axial size of the atomic ensemble 
(~ 2 mm FWHM), the spatial resolution of the in-situ 
absorption imaging with the CCD camera was not high 
enough mostly because of the aberrations of the imaging 
lens to resolve the radial size, resulting in an underesti¬ 
mation of the atom density. 

In a later independent experiment, hereafter referred 
to as a ‘reference experiment’, we estimated the radial 
size from the temperature of the atomic ensemble and 
the trap frequency of the radial direction, setting the 
trapping conditions almost the same as those employed 
in the Ramsey measurements. Those trapping conditions 
set almost the same were (1) the same loading sequence 
as described in the preceding subsection; (2) the dipole- 
trap laser power (with an error described in the next 
sentence); (3) the dipole-trap laser focusing; and (4) the 
total number of atoms (with an error described in the 
next sentence). We performed two reference experiments 
under two different trapping conditions corresponding to 
the higher and lower densities in the Ramsey measure¬ 
ments, setting the dipole-trap laser power (the trapping 
condition (2)) with a difference of ~ 2 % and ~ 5 % from 
the power averaged over the Ramsey measurements for 
the higher and lower densities, respectively, and setting 
the total number of atoms (the trapping condition (4)) 
within 1 s.d. ~ 14 % and ~ 33 % of the eight and seven 
corresponding values measured in the Ramsey measure¬ 
ments for the higher and lower densities, respectively. In 
these reference measurements, the temperature was mea¬ 
sured by an expansion of the atomic ensemble and para¬ 
metric heating was employed to infer the radial trap fre¬ 
quency. The temperatures and the trap frequencies were 
thus obtained to be ~ 67 /iK and 2.2 kHz and ~ 39 /iK 
and 1.1 kHz under those two trapping conditions, respec¬ 
tively. These temperatures and trap frequencies gave 
typical radial sizes (FWHM) of the atomic ensembles in 
the reference experiments to be ~ 14 and 20 /rm, which 
we regarded to be the radial sizes of the atomic ensembles 


in the Ramsey measurements for the higher and lower 
densities, respectively. 

The total number of atoms was also slightly underesti¬ 
mated by the in-situ absorption imaging in the Ramsey 
measurements because of the spatial resolution, so that 
it was calibrated in a later independent experiment in 
which the loading sequence and the trap laser focusing 
(the trapping conditions (1) and (3)) were almost the 
same as those employed in the Ramsey measurements. 
In this calibration experiment, we measured the total 
number of atoms by the absorption imaging with the 
expansion of the atomic ensemble as a function of the 
one measured without the expansion to obtain a linear 
calibration curve with its slope being 1.04 ± 0.03. This 
linear calibration curve gave the total numbers of atoms 
to be ~ 6 X 10® and ~ 4 x 10^ for the higher and lower 
densities in the Ramsey measurements, respectively, for 
n = 42. 

These radial sizes and the total numbers of atoms were 
combined with the axial sizes measured in situ in the 
Ramsey measurements, to give the peak atom densities 
of ~ 1.3 X 10^^ and ^ 4 x 10^® cm“® for the higher and 
lower densities, respectively, for v = 42. Similarly, the 
peak atom density and the total number of atoms for 

= 38 were estimated to be ~ 1.2 x 10^^ cm“® and 
~ 5 X 10® for the higher density and ~ 4 x 10^° cm“® 
and ~ 4 X 10"^ for the lower density, respectively, and for 

= 50 they were estimated to be ~ 1.2 x 10^^ cm“® and 
~ 4 X 10® for the higher density and ~ 3 x 10^® cm“® and 
~ 3 X 10"^ for the lower density, respectively. 


Estimation of the temperature 

The temperature of the atomic ensemble was not mea¬ 
sured in situ in a series of the Ramsey measurements; 
however, it was measured in a later independent exper¬ 
iment, hereafter referred to as a ‘temperature experi¬ 
ment’, in which we set the trapping conditions (1)(4) 
in the same way as described in the preceding subsec¬ 
tion. In this temperature experiment, the temperature 
was measured by an expansion of the atomic ensemble. 
Two temperature experiments were performed indepen¬ 
dently under the trapping conditions corresponding to 
the higher density in the Ramsey measurements, giving 
~ 67 and ~ 72 /xK, respectively, so that we estimated 
the temperature of the higher density ensemble to be 
~ 70 /iK. Similarly, a temperature experiment was per¬ 
formed independently under the trapping conditions cor¬ 
responding to the lower density in the Ramsey measure¬ 
ments, giving ~ 39 /iK, so that we estimated the temper¬ 
ature of the lower density ensemble to be ~ 39 /tK. 
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Rydberg excitation and detection 

The output of a Tirsapphire laser system (Spectra 
Physics; Spitfire Ace, wavelength ^ 779 nm, pulse width 
~ 1 ps, repetition rate 1 kHz) was used as the infrared 
pulse and was also used to pump an optical paramet¬ 
ric ampliher (Spectra Physics; TOPAS) to generate the 
blue pulse tuned to ^ 481 nm. The repetition rate was 
reduced using pulse pickers to synchronize it with the 
atom preparation sequence whose duration was 1.6 s. To 
reduce the number of Rydberg states to be excited, the 
spectra of those infrared and blue pulses were cut using 
homemade pulse shapers in a 4/-setup (/ = 500 mm), 
respectively. The infrared and blue pulses were com¬ 
bined collinearly with a dichroic mirror and their rela¬ 
tive timing was coarsely adjusted to be zero by cross¬ 
correlation measurements based on sum frequency gen¬ 
eration and was further optimized to maximize the Ryd¬ 
berg ion signals. They were introduced into a Michelson- 
type interferometer to produce a pair of identical dou¬ 
ble pulses, each of which was composed of the infrared 
and blue pulses. Those two double pulses induced the 
pump and probe excitations, respectively (see Fig. 

The relative phase between these two double pulses were 
tuned with attosecond precision [IS]. Those two dou¬ 
ble pulses were combined collinearly with the dipole-trap 
laser beam with another dichroic mirror. Those two dou¬ 
ble pulses and the trapping laser beam were focused with 
a plano-convex lens (/ = 250 mm) to the atomic ensem¬ 
ble. 

The infrared and blue pulses, and the optical pump¬ 
ing beam were circularly polarized in the same direc¬ 
tion with respect to the magnetic field, so that the state 
z/D 5 / 2 ,mj = -1-5/2 was mostly populated, and excitations 
to S Rydberg states were suppressed due to transition 
selection rules even though the effective two-photon ex¬ 
citation spectrum covered the S states. This excitation 
scheme also suppressed the Raman transition between 
the F=2 and F=1 hyperfine states in the ground state 
that could happen within the single infrared pulse to in¬ 
duce undesirable beating in the Ramsey contrast as a 
function of r, of which an example is shown in Fig. ii, 
with a period of 146 ps, which is the reciprocal of the 
hyperfine splitting 6.83 GHz. 

The maximum population of the lyD Rydberg states 
was not > 5 %, to suppress photo-ionization. The typical 
pulse energies of the infrared and blue pulses were ^10 
and 400 nJ, respectively, for the 1.2% population and 
were ^ 30 and 600 nJ, respectively, for the 3.3% popula¬ 
tion. The population was estimated from the loss of the 
number of the ground-state atoms in the higher-density 
ensemble induced by the irradiation of the picosecond 
infrared and blue pulses. We measured the number of 
the atoms by absorption imaging with and without the 
picosecond pulses and compared those numbers to eval¬ 


uate the loss. See the previous subsection ‘Estimation of 
the atom density’ for the details of the evaluation of the 
atom number by absorption imaging. The loss induced 
by a single pair of the infrared and blue pulses was too 
small (a few percent) to be evaluated securely, so that 
we shined 30 or 50 pairs of the infrared and blue pulses, 
each of which was accompanied by held ionization m, 
at a 1 kHz repetition rate to induce the loss that was 
clearly visible and was used to infer the loss induced by 
the single pair of the infrared and blue pulses. The loss 
for different amounts of pulse pairs followed the relation 
(A^r/At) = (1 — Pe)'^, where N,- is the number of atoms 
remaining after q pulses and At is the total atom num¬ 
ber. pe is the Rydberg state population for the atoms in 
the ensemble. 

After the probe excitation, the populated Rydberg 
states were ionized by means of held ionization. The Rb’'" 
ions thus produced were detected with a micro channel 
plate placed 5.5 cm away from the atomic ensemble. The 
electric held for the ionization was triggered 50 ns after 
the probe excitation, reaching the ionization threshold 
within the next 100 ns. The output of the micro chan¬ 
nel plate was amplihed with a preampliher and sent to 
a gated integrator, whose output was fed into a com¬ 
puter. It is noteworthy that the held-ionization spectrum 
shown in Fig.[^ was measured using an oscilloscope with 
the electric held ramped up slowly on the microsecond 
timescale to check how many Rydberg states were pop¬ 
ulated by our broadband excitation with the picosecond 
pulses. 

Time-domain Ramsey interferometry 

The interaction among Rydberg atoms was observed 
by time-domain Ramsey interferometry with a pair of 
two-photon excitations: pump and probe. Their delay r 
was coarsely tuned on the picosecond timescale with a 
motorized mechanical stage placed in one arm of the in¬ 
terferometer mentioned above and was scanned hnely on 
the attosecond timescale with a piezoelectric transducer 
to measure Ramsey interferograms. A He-Ne laser beam 
was introduced to the interferometer to check the linear¬ 
ity of the scan by monitoring its optical interference. The 
period of this optical interference was also used for the 
calibration of the pump-probe delay t. We measured the 
held-ionization signals of the two atomic ensembles with 
different densities alternately to suppress systematic un¬ 
certainties, scanning r in steps of ~ 30 as over a range 
of ^ 3fs at each coarse delay tuned by the mechanical 
stage on the picosecond timescale. The obtained Ramsey 
interferogram was htted with a sinusoidal function. As 
the expected energy shift induced by the interaction was 
at most on the order of 10 GHz, much smaller than the 
eigenfrequency of the Rydberg state itself (^ lx 10^® Hz), 
we htted the interferograms with the same eigenfrequen- 
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cies for the higher and lower densities to evaluate the 
phase shift of the higher-density ensemble from the lower- 
density one. We have defined the contrast of the inter- 
ferogram to be the ratio of the amplitude of the fitted 
sinusoidal function to its mean value. In Fig. |^-c, the 
signal intensities are normalized by the mean value of the 
sinusoidal function fitted to each interferogram. We have 
defined Ramsey contrast to be the ratio of the contrast of 
the higher-density ensemble to that of the lower-density 
one, as shown in Figs and We assume that the 
interactions are negligible in the lower-density ensemble, 
which is thus taken to be a reference for measuring the 
contrast decay and the phase shift that may be induced 
by the interactions in the higher-density ensemble. 

Outline of the exactly solvable model simulation 

We represent each atom as a two-level system, which 
is a pseudo-spin system, consisting of a ground state 
|g) (= 14,)) and an excited Rydberg state |e) (= |t)) 
with energies £'g and Eg, respectively. This assumption 
should be reasonable, as the contribution of the neigh¬ 
bouring Rydberg levels is small for the 42 D 5/2 state, as 
seen from the field-ionization spectrum shown in Fig. [^. 
The Watom wave function |4'(t))jv is initially assumed 
to be a product of independent single-atom wave func¬ 
tions |'I'(0))jv = |g)®'^- The experiment consists of four 
stages: (i) the pump excitation, (ii) an evolution with 
an A^-atom Hamiltonian H, (iii) the probe excitation 
and (iv) the population measurement by field ioniza¬ 
tion. The experimental observable is the number of Ry¬ 
dberg excited atoms detected as a time-dependent sig¬ 
nal P{t) = ('I'(T)|jv.P|'I'(r)) 7 v, which is the expectation 
value of the sum of projection operators P = 
where Pj = |t)j (tij /-^ measures the population in the 
Rydberg state for atom j normalized to the total atom 
number N. 

The iV-atom Hamiltonian in stage (ii) includes the 
atomic energies and the interactions U{rj]^) between 
atoms in the Rydberg states. This is justified on the 
picosecond timescale of the current measurement, as Ry¬ 
dberg interactions provide coupling strengths on the or¬ 
der of GHz for micrometre separations, much larger than 
those between Rydberg- and ground-state atoms and be¬ 
tween two ground-state atoms on the order of kilohertz 
and less than kilohertz, respectively, for micrometre sep¬ 
arations. This leads to a diagonal Hamiltonian of the 
form in equation ([^. 

Continuum approximation 

In this approximation we assume the density distribu¬ 
tion to be homogeneous in a small volume around a par¬ 
ticular position r (ref. 54). The Ramsey signal P{y,t) is 


calculated for such a homogeneous region as 

P(r,r) « 2pgPe3?{l + e*(-"+«(pg+Pe 7 (^))'^“('')-'}, 

( 6 ) 

where 7 (t) = 3^3 drPe'^~^'^ and rn is referred to 

as a cutoff radius, so that we consider interactions only 
among the atoms within a sphere whose radius is tq. 
A blockade radius re is determined by the bandwidth 
of the pump excitation and is ~ 0.88 /rm in the present 
experiment. iVo(r) is the number of atoms in the local 
volume ^(^Q —Tg) and depends on the local density n(r) 
due to A^o(r) = ^(?'o ~ This signal is averaged 

over the whole ensemble whose atom-density distribution 
is assumed to be Gaussian with the peak atom-density 
Up to obtain equation in which the Ramsey contrast 
and phase shift can be derived from 

-y^TTTT-p jQ y \ n / 

(7) 

to be the absolute value and phase of the function g(r), 
respectively. In the case that the potential is given by a 
van der Waals interaction and in the limit Ngin) —> 00 , 
the continuum approximation leads to the expression in 
equation ([^ for the function ^(r). We expect this ap¬ 
proximation to describe the system well for a sufhciently 
low population of Rydberg excited atoms Pe ^ Pg (ref. 
54), which is safely satisfied in our experiments. 
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Supplementary Fig. 1: Two-atom dipolar potentials around the 42 D 5/2 + 42 D 5/2 asymptote, — 0 and M = 

Faji + mj 2 = 5, where 6 is the same angle as has been introduced in Supplementary Eq. (19l. Only the asymptotic states 
with electronic angular-momenta below I = 5, principal quantum-numbers ranging from 38 to 46, and within ±30 GHz from 
the 42D5/2 ± 42D5/2 asymptote are considered in the diagonalization of the potential-energy matrix, (a) The potentials 
that dominate the interaction between two I42D5/2, mj = 5/2) atoms are indicated by black solid lines, (b) The color 
code indicates the fraction of the population of the asymptotic I42D5/2, mj = 5/2) (g) I42D5/2, mj = 5/2) state that is 
contained in an interaction-induced mixed state as a function of the interatomic distance. Only the contribution of the I42D5/2, 
mj = 5/2) g) I42D5/2, mj = 5/2) asymptotic state needs to be considered since it is the only state in the displayed energy range 
that can be addressed by the pulse excitation due to selection rules. 
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Supplementary Fig. 2: Mean-field analysis of the Ramsey contrast and phase-shift with an anisotropic van 
der Waals potential. The black-diamond-shaped and blue-circle data-points show the Ramsey contrasts (a and b) and 
the phase-shifts (c and d) measured with the population of the 42D5/2 state being ~1.2% (a and c) and ~3.3% (b and 
d), respectively. In a and b, the Ramsey contrasts simulated by the mean-field model with the anisotropic potential given 
by Supplementary Eq. ( |19[ ) (black and blue solid lines) are compared with the measured ones. Similarly the measured and 
simulated phase-shifts are compared in c and d. The interaction strength is limited below 75 GHz, which is the bandwidth 
(half width half maximum) of the pump excitation, and the peak atom density is set to ~ 1.3 x 10^^ cm“® in these simulations. 
The coefficient Ce = 3.4 GHz ^m® has been used in these simulations. The error bars represent the standard deviation. 
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Supplementary Fig. 3: The mean-fleld analysis of the Ramsey contrast and phase-shift with the dipole-dipole 
(DD) interaction and the hybrid form of a dipole-dipole and a van der Waals (DD-vdW) interaction without 
anisotropies. The black-diamond-shaped and blue-circle data-points show the Ramsey contrasts (a and b) and the phase- 
shifts (c and d) measured with the population of the 42 D 5/2 state being ~1.2% (a and c) and ~3.3% (b and d), respectively. 
In a and b, the Ramsey contrasts are simulated by the mean-field model with the DD (dahed line) and DD-vdW (dotted line) 
interactions without anisotropies and are compared with the measured ones. Similarly the measured and simnlated phase-shifts 
are compared in c and d. The mean-field simulations with the pure van der Waals (vdW) interaction (solid line), which have 
been shown in Fig. 3 in the main text, are presented again to be compared with the DD and DD-vdW result. The interaction 
strength is limited below 75 GHz, which is the bandwidth (half width half maximum) of the pump excitation, and the peak atom 
density is set to ~ 1.3 x 10^^ cm“® in these simulations. The coefficient C 3 = 1 GHz /rm® has been used for the DD interaction, 
and a combination of C 3 = 3.4 GHz ^m^ and = 0.81pm have been used for the DD-vdW interaction [see Supplementary 
Eq. (20l]. The error bars represent the standard deviation. 
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Time delay from the pump excitation (ps) 

Supplementary Fig. 4: State-resolved field-ionization spectra, (a) Ion signal measured as a function of time-delay from 
the ps pnlsed-laser excitations to Rydberg states v = 38. The electric field was ramped at the same speed as the one for Fig. Ic 
in the main text. The atom density and Rydberg population were ~ 4 x 10^° cm~^ and 3.2 ± 0.1 %, respectively, (b) The ion 
signal for = 42 with the atom density and Rydberg popnlation being ~ 4 x 10^° cm“^ and 1.2 ± 0.1%, respectively. This 
is the same signal as shown in Fig. Ic. The red-shaded region indicates the integration range nsed for the estimation of the 
relative population in the 43D state, whereas the blue-shaded area shows the integration range for the 42D and 41D states, (c) 
The ion signal for u = 50 with the atom density and Rydberg population being ~ 3 x 10^° cm“^ and 3.1 ± 0.2 %, respectively. 
The ramp-up speed of the electric field was slower than that for v = 42 and 38 by a factor of 2/3. The bandwidths of the 
excitations in these measurements in a - c are the same as the one for the pump excitation in the Ramsey measurements. 
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Supplementary Fig. 5: Convergence of the simulated Ramsey contrast and phase-shift as functions of an average 
number of interacting atoms. The Ramsey contrast (a) and phase-shift (b) at r = 500 ps are simulated by the theory 
model with the continuum approximation and are plotted as functions of the cutoff radius ro (the lower abscissa) and of an 
average number of interacting atoms within the volume V = ^(I'o ~ i"b) (the higher abscissa). The population of the 42 D 5/2 
is set to ~ 3.3 % in these simulations. The interaction strength is limited below 75 GHz, which is the half width half maximum 
of the pump excitation, and the peak atom density is set to ~ 1.3 x 10^^ cm“® in these simulations. The red dashed and green 
dotted lines show the results with the dipole-dipole interaction and the hybrid form of a dipole-dipole and a van der Waals 
interaction without anisotropies, respectively. The results with the van der Waals interaction without an anisotropy, which has 
been used in the main text, are displayed by the blue solid lines. The dark-grey solid lines represent the measured Ramsey 
contrast and the phase-shift, each of which is the average over eight points around r = 500 ps in Figs. 6b or d. The light-grey 
shaded area represents one standard deviation of the average over those eight measured values. The black solid line shows 
the Ramsey contrast |g(r)| given by Eq. (7) in the main text with 7 (r) = 0, giving the upper limit of the contrast decay and 
accordingly the lower limit of the number of atoms ~ 32 to reproduce the Ramsey contrast ~ 0.45 measured at t = 500 ps, 
irrespective of the potential curves. 
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Supplementary Fig. 6 : Numerical simulation of the photo-ions contribution to the Ramsey signal, (a) Ramsey- 
contrasts simulated as functions of pump probe delay r in the presence of photo-ions for u = 38 (black), 42 (blue), and 50 (red), 
(b) Phase-shifts simulated as functions of pump probe delay r in the presence of photo-ions for v = 38 (black), 42 (blue), and 
50 (red). 




Supplementary Fig. 7: Simulation of the effects of the Rydberg interactions during the excitation on the 
Ramsey signal, (a) The Rydberg populations at the end of the excitation pulse as functions of the pulse duration with 
(red) and without (blue) the interaction, (b) the Ramsey oscillations around the pump-probe delay r ~ 400 ps with (red) and 
without (blue) the interaction during the excitation pulses. 
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Supplementary Fig. 8 : High resolution measurement of the Ramsey-oscillation contrasts, (a) The contrasts as 
functions of the pump-probe delay r for the higher-density (red dots) and lower-density (blue dots) ensembles, respectively. 
The red- and blue-shaded parts are numerically simulated recurrence-motions of the Rydberg wave-packet for the higher- and 
lower-density ensembles, respectively, (b) Field ionization spectrum associated with the measurement shown in a. 
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Supplementary Fig. 9: The calibration of the peak atom density as a function of the power of the dipole-trap- 
laser. The red-circle data points show the peak atom density measured in the calibration experiment as a function of the 
dipole-trap-laser power. The blue-diamond-shaped data and their error bars are not measured directly, but are obtained by 
linear interpolation. These blue-diamond-shaped data points represent the estimated peak atom densities used in the density 
dependence measurement of the Ramsey contrast shown in Fig. 4b in the main text. The error bars represent one standard 
deviation. 
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Supplementary Fig. 10: The principal-quantum-number dependences of the phase-shift, (a - c) Measured phase- 
shifts are plotted as functions of r for three different Rydberg levels u = 38, 42, and 50. Each of them has been measured 
simultaneously with the Ramsey contrast with the same principal-quantum-number shown in Fig. 4a in the main text. The 
simulations indicated by the black, blue, and red solid-lines have been performed by the theory model with the continuum 
approximation for the van der Waals interaction with the adjusting parameters being Cq — 8GHz/^m® for u = 38, Cq — 
34GHz/im® for u = 42, and Cq = 103 GHz ^m® for u = 50. The peak atom-density is set to the estimated density for each 
Rydberg level in these simulations. The density estimations are described in Methods section “Estimation of the atom density” 
in the main text. It should be noted that several Rydberg states are excited in the case of u = 50, whereas we have considered 
an excitation only to a single Rydberg state to perform the simulations for all of the three Rydberg levels. The error bars 
represent the standard deviation. 
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Supplementary Fig. 11: The theory-model analysis of the Ramsey contrast and phase-shift beyond mean-field 
with the continuum approximations and with an anisotropic van der Waals potential. The black-diamond-shaped 
and blue-circle data-points show the Ramsey contrasts (a and b) and the phase-shifts (c and d) measured with the population 
of the 42 D 5/2 state being ~1.2% (a and c) and ~3.3% (b and d), respectively. In a and b, the Ramsey contrasts simulated with 
the anisotropic potential given by Supplementary Eq. ( |19| l (black and blue solid lines) are compared with the measured ones. 
Similarly, the measured and simulated phase-shifts are compared in c and d. These simulated results have been obtained by 
Eq. (6) in the main text combined with the potential anisotropy, employing the cutoff radius ro ~ 4 /rm and its corresponding 
atom-number Nq = 450 at the peak density. The interaction strength in these simulations is limited below 75 GHz, which is the 
half width half maximum of the pump excitation, and the peak atom density is set to ~ 1.3 x cm“® in these simulations. 
The coefficient = 63 GHz /rm® has been used in these simulations. The error bars represent the standard deviation. 
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Supplementary Fig. 12: The theory-model analysis of the Ramsey contrast and phase-shift beyond mean-field 
with the continuum approximation and with the dipole-dipole (DD) interaction and the hybrid form of a 
dipole-dipole and a van der Waals (DD-vdW) interaction without anisotropies. The black-diamond-shaped and 
blue-circle data-points show the Ramsey contrasts (a and b) and the phase-shifts (c and d) measured with the population 
of the 42 D 5/2 state being ~1.2% (a and c) and ~3.3% (b and d), respectively. In a and b, the Ramsey contrasts simulated 
with the DD (dahed line) and DD-vdW (dotted line) interactions without anisotropies are compared with the measured ones. 
Similarly the measured and simulated phase-shifts are compared in c and d. The simulated results for the DD interaction have 
been obtained by Eq. (6) in the main text, employing the cutoff radius ro ~ 4 /rm and its corresponding atom-number No = 450 
at the peak density. It should be noted that these DD results are pushed down as the cutoff radius is further increased, as is 
indicated by arrows in the hgures and will be further discussed with Supplementary Fig. The simulations with the pure van 
der Waals (vdW) interaction (solid line), which have been shown in Fig. 6 in the main text, are presented again to be compared 
with the DD and DD-vdW results. The interaction strength is limited below 75 GHz, which is the bandwidth (half width half 
maximum) of the pump excitation, and the peak atom density is set to ~ 1.3 x 10^^ cm“® in these simulations. The coefficient 
C 3 = 2.5 GHz ^m^ have been used for the DD interaction, and a combination of C 3 = 4.3 GHz ^va? and = 1.93 /rm have been 
used for the DD-vdW interaction [see Supplementary Eq. (20l]. The error bars represent the standard deviation. 
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Supplementary Fig. 13: Phase-shifts calculated as functions of the pump-probe delay r. (a) Comparison between the 
exact (solid lines) and mean-field (dotted lines) solutions with the Gaussian atom-density distribution for different populations 
Pe (~ 1-2, ~ 3.3, and ~ 10.0%) in the 42D Rydberg state. The Gaussian distribution has been taken from the Ramsey 
experiment shown in the main text. The Ce coefficients are set to 34GHz/rm® for the exact solution and 1.9 GHz pm® for the 
mean-field solution, respectively, which are the values used in the calculations performed in the main text, (b) Comparison 
between phase-shifts calculated with the Gaussian (solid lines) and homogenous (dashed lines) atom-density distributions within 
the mean-field model for different populations Pe ~ 1.2% and ~ 3.3% in the 42D Rydberg state. The Gaussian distribution is 
the same as in a, whereas the homogeneous density has been set to the average density of the Gaussian distribution. The Ce 
value is set the same as in a. 


Supplementary Notes 

Supplementary Note 1: Estimation of the photo-ion contribution to the Ramsey signal 

In order to investigate the effect of the electric field generated by photo-ions on the Ramsey signal, we have 
performed numerical simulations for v = 38, 42, and 50, including ion contributions. In each of those simulations, 
atom locations in an ensemble with a Gaussian distribution are generated by the Monte-Carlo method. In this list of 
atom locations a certain number of random locations are taken to be ion locations. The number of ions is obtained as 
follows. The ion fractions are estimated for zz = 38 and 50 from the field ionization spectra shown in Supplementary 
Figs. & and c, respectively, and for zz = 42 from the spectrum similar to the one shown in Supplementary Fig. to 
be the ratios of the areas of the ion peaks at ~ 5.6 p,s for i/ = 38, 42 and ~ 6.2 /is for zz = 50 to the total areas of the 
spectra. The numbers of ions are thus estimated as the products of those ion fractions, our Rydberg populations Pe’s, 
and the total number of atoms to be 109 for zz = 38 with pe ~ 3.2%, 232 for zz = 42 with pe ~ 3.3%, and 108 for 
z^ = 50 with Pe ^ 3.1 %, respectively. These numbers are the upper limits of the number of photo-ions as they may 
include the ions produced at much later time (100 ns timescale) by Rydberg interactions [1]. At each atom position 
in each of the ensembles with these ion numbers, the joint electric field generated by those ions is calculated. The 
Stark shift induced by this joint electric field is calculated for each atom, so that the period of its Ramsey oscillation 
is shifted accordingly. Those Ramsey oscillations are averaged over the whole ensemble of the Rydberg atoms to give 
the contrast-decay and phase-shift shown in Supplementary Fig. |^for each oi u = 38, 42, and 50. It is seen that 
the contrast-decay and phase-shift are less than 2% and 2.5 degrees at the longest pump-probe delay 500 ps in our 
measurements and are almost negligible. 

Supplementary Note 2: Discussion on the pump-probe excitation and the zero-delay 
offset in the phase-shift 

We consider the time-domain Ramsey interferometry with the identical pump and probe excitations, each of which is 
a one-photon excitation of a two-level atom, consisting of a ground state |g) and an excited state |e), with rectangular 
laser pulses for simplicity. The results can be easily adapted to near-resonant two-photon excitations that have been 
used in the present experiments. 
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The time dependent Schrodinger equation is given by 

= h\'4)) 

= (^||e)(e| - ^||g)(g| - deg • E(t)|e)(g| - d*g • E(t)|g)(e|) |^/>). 

Here, w is the angular frequency for the atomic transition, E(t) = Eo(0(t — to) — 0(t — ti)) cos(a;;(t — to)) is the 
pump laser field where ui\ is the laser frequency, 0(t) is the Heaviside step function, to is the beginning of the pulse 
excitation, ti is the end of the excitation, \ijj) = C'g(t)|g) + C'e(t)|e) is the atomic state, and dog = (e|d|g) is the 
dipole-matrix element. Solving the Schrodinger equation Eq. Q by employing the rotating-wave approximation for 
the evolution of the atomic state during the pump excitation we obtain 


/Cg(ti)\ ^ / (cos(6»/2)-ir?sin(6l/2))e*T<5‘ -i^sm{e/2))e^^e^^^* \ C'g(io) \ 

V C'e(H) / y —i^sin( 0 / 2 )e“*^e“®" 2 *'^* (cos(0/2)-|-i?; sin(0/2))e“*'^^‘y V C'e(io) / 



with 6 — flSt, 77 = A/0, ^ = Or/O, and O = y^Op^ + A^, where Or is the Rabi frequency, and A = w — wi is the 
detuning of the laser frequency uj\ from the atomic resonance lo. The time 5t = ti — tg is the duration of the pump 
excitation, and y represents the phase of the complex Rabi frequency Or = Orc*^. It is important to note in the 
formulation given in Eq. 0 that the matrix H, which describes the time evolution of the atomic state during the 
pump excitation, is given in the lab frame, in which the field-free evolution of the atomic state between the pump 
and probe excitations is described. 

The matrix for the excitation can be expressed in the form 




Icgle®"^/^ 7|ce|e*^/^ \ 


( 2 ) 


where (j) and C are phases of the complex amplitudes Cg and Ce, respectively. The pump and probe excitation operation 
on the ground state |g) results in 


|i/)(t)) = Hexp(-jiJr/;i)H|g) (3) 

|cg|V^e*t--|c3|2e-t- \ 

*|cg||ce|e-*«/2(e*'^/2e*t^ -F J 

Therefore, the population in the excited state P{t) = ('(/;(T)|e)(e|'(/(T)) that remains after the probe excitation is given 
by 



P{t) = 2|cg|^|cep (1-I-cos(a;r-I-/))). (4) 

Here only the phase (j) remains in the Ramsey signal, while the phase ( is cancelled out. According to Eq. Q and 
Eq. § the phase is defined by 

r2sin(wi(5t/2) cos(ni5t/2) — Acos(a;ii5t/2) sin(fI(5t/2) 

H cos(a;i5t/2) cos{il6t/2) + A sin(a;ii5t/2) sin(r2(5t/2)' 


When the laser detuning vanishes (A = 0) we obtain = uJiSt = LoSt. Similarly, if Hr —>• 0, the phase (j) converges 
to uiSt. In the regime with |A| >> |Hr| where H = y -|- A^ « A -|- the effects of the laser field can be 

q 2 

characterized by the AC-Stark shift Wac = yj, so that we obtain 


tan(^/2) 



Therefore the phase aquired during the excitation is given by the AC-Stark-shifted phase-evolution 


( 6 ) 


(j)Ki {cO - 20Ja,c) St. 


( 7 ) 
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If the interactions among the atoms are neglected during the ps pulses, this intrapulse phase is common to the higher- 
and lower-density ensembles under the condition that their sizes, shapes, and positions are identical. Accordingly, the 
intrapulse phase is cancelled out in the phase-shift between those two ensembles. 

In the actual measurements, however, slight differences between the sizes, shapes, and positions of the higher- and 
lower-density ensembles can lead to different intrapulse phases due to different AC-Stark shifts (see Methods section 
“Estimation of the atom density” in the main text). Additional intrapulse phases due to the Rydberg interactions 
during the pump and probe excitations are considered to be negligibly small as seen in the theoretical simulations 
demonstrated in the next paragraph. 

We have performed theoretical simulations of the Ramsey signals to investigate the effects of the Rydberg in¬ 
teractions during the pump (or probe) excitation, whose width is 10 ps, on the contrast and phase of the Ramsey 
oscillations. We consider a two-level system with its level spacing that corresponds to a frequency of 1.008 x 10^^ Hz 
and is the energy difference between 5S and 42D states. We assume the van der Waals interaction Cg ~ 34 GHz /im® 
as it has been obtained in the beyond-mean-field analyses of our data in the main text. The interatomic distance 
that gives the interaction timescale comparable to the width of the excitation ^ 10 ps is estimated to be ^ 0.84/rm. 
The average number of atoms within a sphere whose radius is 0.84 /im is estimated to be less than two in our higher- 
density ensemble. The number of atoms N that has been considered in the present simulation is five and is thus large 
enough to investigate the interaction effects within the ~ 10 ps pulse excitation. We have compared the populations 
of the excited state generated by an effective off-resonant two-photon excitation with and without the atom-atom 
interactions during the excitation process. We have considered the atoms to be randomly distributed within a cube 
with a volume rg = A^/riav where the average atom density riav = (1.3/2^^/^^) x 10^^ cm^ is the same value as the one 
in our experiment. We plot the excited state populations at the end of the excitation pulse as functions of the pulse 
duration with and without the atom-atom interactions in Supplementary Fig. Here the red and blue curves show 
the results with and without the interactions, respectively, during the excitation. We fixed the average Rabi frequency 
to 18.88 GHz, and the detuning from resonance emerging due to the two-photon light shift to 150 GHz, so that the 
resulting populations are ^ 3 % in the excited state at the pulse duration of 10 ps. It is seen in Supplementary Fig. [Tfr 
that the populations are almost the same between these two cases with their population difference being less than 
1 % at 10 ps, indicating that the interaction effects on the excited-state population are almost negligible during the 
pulse excitation. Next we have simulated the Ramsey signals with the pump and probe excitations and evaluated 
the intrapulse-interaction effects on the contrast and phase of the Ramsey oscillations. The simulated Ramsey os¬ 
cillations are shown in Supplementary Fig. where the red and blue curves show the results with and without the 
interactions during the excitation pulses, respectively, around the pump-probe delays r = 400 ps. We calculate the 
contrasts for those traces and obtain 0.867 and 0.844 with and without the interactions, respectively. The difference 
between the contrasts is thus ^ 3 %, originating from the interaction during the the excitations. The phase difference 
between those two traces is less than ~ 1 degrees. This corrsponds to the intrapulse phase discussed in the preceding 
paragraph. These simulated values for the contrast reduction and the phase difference are much smaller that the 
measured contrast-decay and the amount of phase-shift seen in Figs. 3 and 6. It is thus concluded that the interaction 
effects on the contrast and phase are almost negligible. It is also concluded that the contribution of the interactions 
to the intrapulse delay is negligibly small. 

We have verified that the discussion above holds entirely also for a three-level system with rectangular laser pulses. 
We consider a system consisting of a ground state |g), an intermediate state |i), and a Rydberg excited state |e). 
A laser field with frequency wn couples the state |g) to the upper state |i) with Rabi frequency and detuning 
Ai = wii — Wi, where Wi is the corresponding transition frequency. A second laser field with Rabi frequency drives 
the transition between the states |i) and |e) with the transition frequency uj 2 - In the case that |Ai| ^ |Hpji|, |Hn 2 |) we 
can neglect the excitation to the intermediate state, so that the three-level system is effectively reduced to a two-level 
system. In our actual experiments, however, we use an excitation laser-pulse whose envelop is not rectangular. The 
condition |Ai| ^ |f^Ri|, |^^R 2 | holds in the beginning and end of the pulse, whereas it might be violated at its peak 
intensity around the middle of the pulse, so that the intermediate state |i) can be populated at the peak intensity. 
However, this state |i), which is not a Rydberg state, has a negligibly small intrapulse-interaction effect during the 
excitation pulse and due to adiabaticity given by the condition wi, {uj 2 — wi) ^ |f^Ri|, |f^R 2 | no population remains 
in the intermediate state |i) after the pulse. This shows that our two-photon excitation can be effectively treated 
as a one-photon excitation. The discussions above on the intrapulse phase and the interaction effects during the 
excitation also apply to the two-photon excitation in our current experiments. 

Supplementary Note 3: Field ionization spectra and estimation of excitation bandwidth 

In order to reduce the number of Rydberg states to be excited, the spectra of the ps IR and blue pulsed lasers were 
cut to be about 0.13 nm and 0.20 nm, respectively, with pulse shapers in a 4f configuration (/ = 500mm). These 
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bandwidths were common to the Ramsey measurements of the three different Rydberg states v = 38, 42, and 50. 
Supplementary Figure shows state-resolved field-ionization spectra measured by ramping the electric field slowly 
on the microsecond timescale. The spectra indicate that a single state was predominantly populated for each of the 
excitations to = 38 and 42 as seen in Supplementary Figs. and b, whereas more levels were populated for the 
excitation to zz = 50 as seen in Supplementary Fig. This is because the energy levels of Rydberg states are inversely 
dependent on and more congested for higher states. Due to the characteristics of field ionization, the threshold 
value of the electric field to induce ionization depends on v~^ [2], so that higher states are less resolvable with the 
same ramp-up speed of the electric field. Therefore the ramp-up speed was slower for z/ = 50 than for z/ = 38 and 42 
by a factor of 2/3. 

The bandwidth of our Rydberg excitation with the IR and blue pulses was determined from the field-ionization 
spectrum for z/ = 42 presented in Fig. Ic as well as in Supplementary Figs. i). Assuming a symmetrical population 
distribution with respect to the center level z^ = 42, we estimated the relative populations in the excited Rydberg 
states to be ~ 84% for the center state 42D and ~ 8% for each neighboring state 43D and 41D, respectively, from 
the integrated areas of the ion signal peaks. In Supplementary Figs. the integration range for the 43D state (the 
42D and 41D states) is indicated by the red-shaded (blue-shaded) region. We did not deconvolute the ion spectrum 
for those integrations although those peaks were not fully resolved. Based on the relative populations and assuming 
a Gaussian excitation spectrum, we obtained a bandwidth of ~150GHz (FWHM). 

Supplementary Note 4: Structures in the pump-probe delay dependence of the Ramsey 
signal 

We have made another set of supplementary Ramsey measurements with smaller steps of the pump-probe delay r 
than those of the measurements shown in the main text. The results of those supplementary measurements are shown 
in Supplementary Fig. in which the red and blue dots show the contrasts of the Ramsey oscillation for the higher- 
and lower-density ensembles, respectively, as functions of r. It should be noted here that Supplementary Fig. Ep shows 
a field ionization spectrum associated with the Ramsey measurement shown in Supplementary Fig. |^, indicating the 
41D, 42D, and 43D states predominantly populated, similar to the measurements shown in the main text, but more 
fractions of the 41D and 43D states included in these supplementary measurements with a broader bandwidth of the 
ps excitation laser pulse. It is thus reasonable that the recurrence motion of the Rydberg wave-packet with a period of 
^ 10 ps is more pronounced and better resolved in Supplementary Fig. than in Fig. 2d in the main text due to the 
more fractions of the neighboring levels and smaller steps of t, respectively. The red- and blue-shaded parts shown in 
Supplementary Fig. correspond to the recurrence motions of the Rydberg wave-packet numerically simulated only 
with those three Rydberg levels and with the decay factor exp(-a-\/r) seen in Eq. (5) in the main text. It is seen from 
this figure that the simulated results show excellent agreements with the measured r dependences of the contrast, 
demonstrating that our Ramsey signal is not affected by other angular momentum states such as S, P, and F states 
or by the oscillation resulting from the ground-state hyperfine splitting. 

It is understood from the comparison between the structures seen at r ^ 130 — 170 ps in Fig. 2d in the main text 
and Supplementary Fig. that those oscillations can be assigned to the recurrence motion of the wave-packet. It 
is also understood from the comparison between the collapse and revival of the wave packet seen in Supplementary 
Fig. §1, which is due to an anharmonicity of the Rydberg levels, and the global structure seen in Fig. 2d that the 
structures on the ~ 100 ps timescale seen in Fig. 2d are not always assigned to the collapse and revival, but also to 
experimental fluctuations. 

Supplementary Note 5: Active control of many-body dynamics 

Supplementary Note 5-1: Atom-density dependence of the Ramsey-contrast decay 

Figure 4b in the main text shows the Ramsey contrasts for z/ = 42 with a population of 3.5 ± 0.3% at r = 300 
and 510 ps for several different atom densities ranging from the lower to the higher densities described in the main 
text. These measurements were made by changing the power of the dipole-trap laser and thereby the trap depth. 
At first the atom densities in these Ramsey measurements were estimated solely from the total number of atoms 
and the size of the atomic ensemble obtained by in-situ absorption imaging with a GGD camera without expanding 
the atomic ensemble, but were underestimated because of the spatial resolution, as is described in Methods section 
“Estimation of the atom density” in the main text. Therefore it was calibrated in a later independent experiment 
in which the trapping conditions (1) and (3) (the loading sequence and the dipole-trap laser focusing), which were 
described in Methods section “Estimation of the atom density”, were almost the same as those employed in the 
Ramsey measurements. In this calibration experiment, we measured the radial trap-frequency, the temperature, the 
axial size, and the total number of atoms to obtain the atom density as a function of the power of the dipole-trap 
laser. Supplementary Figure shows the results of this calibration measurement accompanied by a calibration curve, 
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which is a linear interpolation function. The atom densities in the Ramsey measurements were estimated from this 
calibration curve as a function of the dipole-trap-laser power as shown by the blue diamond-shaped data points in 
Supplementary Fig. and then the Ramsey contrasts measured at t =300 and 510 ps are plotted against these 
estimated atom densities in Fig. 4b in the main text. It is seen in Fig. 4b that the contrast decay is accelerated as 
the atom density is increased. 


Supplementary Note 5-2: The principal-quantum-number dependences of the phase- 
shift 

Figure 4a in the main text shows the Ramsey contrasts as functions of r for three different Rydberg levels v = 38, 
42, and 50. The populations Pe and estimated peak atom-densities are Pe ~ 3.2 % and ^ 1.2 x 10^^ cm“^ for v = 38, 
Pe ~ 3.3% and ^ 1.3 x 10^^ cm“^ for v = 42, and Pe ~ 3.1% and ~ 1.2 x 10^^ cm“^ for v = 50, respectively (see 
Methods section “Estimation of the atom density” in the main text for these density estimations). It is seen from 
this figure that the dephasing is accelerated by increasing the principal quantum number v of the Rydberg level. We 
have also measured the corresponding phase-shifts for these three levels simultaneously with the Ramsey contrasts 
as shown in Supplementary Fig. |10| The theory-model simulations with the continuum approximation indicated by 
solid lines agree well with the measured results. It should be noted that several Rydberg states are excited in the 
case of ^ = 50, whereas we have considered an excitation only to a single Rydberg state to perform the simulations 
for all of the three Rydberg levels. 


Supplementary Note 6: Calculation of the Ramsey contrast and phase-shift with 
nearest-neighbor interactions 

We follow previous Ramsey studies on Rydberg interactions [3,4] to calculate the Ramsey-contrast-decays and the 
phase-shifts expected for nearest-neighbor interactions. In the Ramsey measurement in the main text, the population 
in the Rydberg state is observed as a function of the delay time r between the pump and probe excitations. Within 
the delay time, nearest neighbor atoms evolve under a Hamiltonian 

2 

H = ^(F;g|g)j(g|j-f F;e|e)j(e|j)+ C(r)|e)i|e) 2 (e| 2 (e|i, (8) 

i=i 

where |g) and |e) are ground- and excited Rydberg-states with energies Eg and Ee, respectively, and U{r) is an 
interaction energy between a pair of Rydberg atoms separated by r. The time-domain Ramsey signal is obtained by 
solving a Schrodinger equation with this Hamiltonian as 


P{t) = 2pgPe (1-l-pgCOs(a;r-I-(())-|-peCOs((w-I-C/(r)//i)T-I-(())) (9) 

= 2pgPe3? 11 -I- ^Pg -I- PeC • 

Here Pg and Pe are the ground- and Rydberg-state populations, respectively, ui = (Ee — Eg)/h is the atomic-resonance 
frequency, and </> is the phase offset arising from the AC-Stark shifts during the pulse excitation. This result is identical 
to the Ramsey signal obtained by Eq. (3) in the main text with N = 2. In a homogeneous atom distribution with a 
density of n, the nearest-neighbor distribution is given by 


Pr{n, r) 


exp 




47rnr^ 


( 10 ) 


where Pr{n,r)dr = 1. The Ramsey signal given by Supplementary Eq. ^ is averaged over this distribution as 
follows: 


poo 

Pav{n,T) = / Pr{n,r)P{T)dr 

Jo 

= 2pgPe5R 11 -h ^Pg -f Pe j 


( 11 ) 


I iU(r)T 

Pr(n,r)dr + / Pr{n,r)e ** dr 
Jru 


where re is the blockade radius determined by the finite bandwidth of the excitation with the IR and blue pulses. 

This is further averaged over the density distribution of our atomic ensemble, which can be modeled by a Gaussian 
distribution to be 

n(x) = ripC 


( 12 ) 
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where rip = ^ 2 ^) 3 / 2 g .2 „ is the peak density, N is the number of atoms in the ensemble, Cxy is the width in the 
x,y-direction, and is the width in the z-direction of the atomic ensemble. The Ramsey signal thus averaged over 
the whole ensemble is given by 


47r 

Rav(r) = — dpp^n{p)P^^{n{p),T) 

dn^\n (^^^Pav(n,T), 


JO 


(13) 


where the radius p is defined by p = ^+ y'^ + ■ In going from the first line to the second line in Supple¬ 

mentary Eq. (13), an integral over the volume of the ensemble is converted to an integral over the density. 

It should be noted that the maximum threshold of the Ramsey-contrast-decay expected for nearest-neighbor 
interactions is determined by the population of a Rydberg state as follows. When the Ramsey oscillations are 
maximally dephased, the third term in the parentheses in the first line of Supplementary Eq. (|^ vanishes. In this 
limit, the oscillation contrast and the phase-shift converge to Pg (= 1 — Pe) and </>, respectively, seen in the second 
term. These features do not depend on the character of the interaction U{r) such as van der Waals and dipole-dipole. 


Supplementary Note 7: Mean-field model 

For the mean-field approximation we follow the notation |e)j(e|j = which results in the following Watom 
Hamiltonian 

N N-1 N 

h = Y. + E E U{n,j)P^^Pg>, (14) 

i=i i=i i>j 

where Eg is the energy of Rydberg state and U{rij) describes the interaction between atoms i and j separated by Vi j. 
Without loss of generality we have set the ground-state energy Eg = 0. By using a relation Pee ^ = (PeP) +SpM^ where 
{Ppp and SpM^ are the mean value and fluctuation of Pp\ respectively, we obtain the mean-field approximation 
from 

P^pP^P « {P^){P^P) + {P^)SP^p + {P^p)SPp^ (15) 

= /p(d\p(I) I /p(j)\p(») _ /p(i)\/p(I)\ 

\-* ee /■* ee ~ \-^ ee /-^ ee \-* ee /\-*^ ee 

which allows to write the Hamiltonian in the following form 

N 

H = (16) 

i=i 

= E ^ + ((^ee) - Ao;,) Pii') - h p(|), 

where w = Eg/H is the atomic-resonance frequency and (pW) is assumed to be common to all the atoms and 
set to (Pee). Each atom can be considered separately, and the interactions enter as shifts of the energy lev¬ 
els with h(^{Pgg) — ^ Aujj and AtUj^ for the Rydberg and ground states, respectively. Here Acjj = 

U{rij)/h is the sum over all interactions with atom ‘j’. 

In this model, therefore, Rydberg interactions modify only the period of the Ramsey oscillation of each atom as 
follows 


Pmf-j{T) = 2pgPg (1 -h COs((w -h PeAuJj)T + (j))) , (17) 

where r is the pump-probe delay, Pg and pg are the ground- and Rydberg-state populations, respectively, (Pee) = Pe, 
and (j) is the phase offset arising from the AC-Stark shifts during the pulse excitation. Oscillations with slightly different 
periods are then averaged over the atom distribution. The interferogram obtained in the Ramsey measurement is 
thus given by 

COs((w -\- PgAljJj)T + (p) 


N 


N 


Pnif{T) — ^ ^ Prai-jiT) — 2pgPg ( 1 + ^ ^ 


i=i 


N 


i=i 


(18) 
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This averaging yields a contrast decay of the Ramsey oscillation. By employing the Monte Carlo method to model 


a realistic distribution of the N atoms as given by Supplementary Eq. (12) we acquire the results of Supplementary 

Eq. pi). 


Supplementary Note 8: Outline of the least-squares fitting 

The least-squares fitting for Fig. 3b in the main text and Supplementary Figs. and was performed so that 
the residual between the Ramsey contrasts measured for ~3.3 % population and their simulations was minimum, 
with reasonable steps of fitting parameters, for 500 atoms, whose configuration was generated by the Monte-Carlo 
simulation. For each of these 500 atoms, the mean-field energy shift is calculated by considering the interactions with 
^ 6 X 10® atoms, and these shifts are averaged over the 500 atoms. The fitting parameters such as the Cq coefficient 
obtained in that least-squares fitting were used to calculate the simulated curves shown in Fig. 3 and Supplementary 
Figs. and with an averaging over 20,000 atoms from the ensemble. 

The least-squares fitting for Fig. 6b in the main text and Supplementary Figs. [TT)d and |12|3 was also performed so 
that the residual between the Ramsey contrasts measured for ~ 3.3 % population and their simulations with Eq (4) 
in the main text was minimum with reasonable steps of fitting parameters. The fitting parameters such as the Cq 
coefficient obtained in that least-squares fitting were used to calculate the simulated curves shown in Fig. 6 and 


Supplementary Figs. 11 and 12 


Supplementary Note 9: Effects of the atom-density distribution on the phase-shifts 

We present detailed analyses of the effects of the atom-density distribution on the phase-shifts. Supplementary 
Figure [13^ shows the exact and mean-field calculations of the phase-shift for the 42D Rydberg state with the Gaussian 
distribution of the atom density. The Gaussian distribution has been taken from the Ramsey experiment shown in 
the main text. Both in the exact and mean-field calculations, the phase-shift is saturated as the Rydberg population 


Pe increases, as shown in Supplementary Fig. 13 r. As the Rydberg population is increased, the initial slope of the 
phase-shift becomes larger due to stronger interactions, so the phase-shift converges earlier to a value around -40 
degrees, and the disagreement between the exact and mean-field results becomes smaller. 

This saturation is due to the Gaussian atom-density distribution of our experimental setup as shown in Sup¬ 


plementary Fig. 13 d, in which the phase-shifts calculated with two different density distributions (Gaussian and 
homogeneous) and two different Rydberg populations are compared within the mean-field model. The Gaussian 
distribution is the same as in Supplementary Fig. [Til, whereas the homogeneous density has been set to the average 
density of the Gaussian distribution. It is seen that the Gaussian distribution gives the phase-shift saturated more 
rapidly for the higher Rydberg population. With the homogeneous distribution, however, the phase-shifts are not 
saturated within this timescale 1000 ps. In contrast to the homogeneous distribution, the rapid decrease of the 
atom density in the Gaussian tails results in the saturation of the phase-shift. This is because the contribution to 
the phase-shift from atoms distant from the center of the Gaussian distribution is suppressed, and therefore the 
phase-shift does not grow afterwards. The stronger interactions for pg ~ 3.3 % yield the phase-shift saturated within 
our measurement time 500 ps. This saturation is reached both by the exact and mean-field calculations for pe ^ 3.3% 
as seen in Supplementary Fig. 13 1 , unlike the case with pe ^ 1.2% and therefore weaker interactions, giving the 
closer agreement between the mean-field and exact results as well as the experimental ones for pe ~ 3.3% than for 
Pe ~ 1.2 %, as seen in Fig. 3d in the main text. 


Supplementary Note 10: Numerical simulations with alternative forms of interactions 
Supplementary Note 10-1: Effective treatment of a two-atom interaction 

At short interatomic distances a dipole-dipole interaction in non-diagonal terms of the Hamiltonian couples the 
initial Rydberg states of a pair of atoms described by | • • • e • • • e • • •) with other Rydberg states | • • • e' • • • e" • • •). 
Such couplings induce hybridization among multiple Rydberg states, leading to congested potential structures. To 
handle this intractable problem in the present study, we have considered a single effective potential for the two-atom 
interaction that represents the congested potential structure. Thereby, the effective interaction enters only the 
diagonal components of the Hamiltonian. 


Supplementary Note 10-2: Numerical simulation with an anisotropic van der Waals 
interaction 

In the main text we have considered an isotropic van der Waals interaction given by U(r) = —C^jr^. Here we 
introduce an anisotropy into the van der Waals interaction by 


C/(r,0) 


Cefl — 3cos^(0))^ 
^6 


(19) 
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where 9 is the angle between the z-axis (Fig. la) and a line connecting two atoms interacting with each other. The 
comparisons between our experimental observations and the simulations with the anisotropic potential U (r, 9) above 
are presented in Supplementary Fig. for the mean-field model and in Supplementary Fig. E] for the theory model 
beyond mean-field with the analytical continuum-approximation, respectively. It is seen in Supplementary Fig. 
that the mean-field simulation fails to reproduce the observed phase-shift again, whereas the theory-model simulations 
agree with the measured Ramsey contrasts and phase-shifts for both of the Rydberg populations Pe ~1.2% to 3.3%. 
The adjustment parameter employed in these simulations with the anisotropic potential is Cq = 63 GHz /rm® and is 
comparable to Cq = 34 GHz /im® employed in the simulations without the potential anisotropy. 


Supplementary Note 10-3: Numerical simulation with a hybrid form of a dipole-dipole 
and a van der Waals interaction and with a pure dipole-dipole interaction 

A hybrid form of a dipole-dipole and a van der Waals (DD-vdW) interaction is defined by 


U{r) 



r < Tc 
r > Tc ’ 


( 20 ) 


where Cc is the crossover radius defined by = Cq/C^. We have performed numerical simulations of the Ramsey 
contrast and phase-shift using this DD-vdW interaction with C 3 and Tc being fitting parameters. 

We have also performed the simulations using a dipole-dipole (DD) interaction without anisotropies. The results 
of those simulations are presented in Supplementary Fig. [^for the mean-field model and in Supplementary Fig. 
for the theory model with the continuum approximation, respectively. The corresponding results with the pure 
van der Waals (vdW) interaction, which have been shown in Figs. 3 and 6 of the main text, are presented again 
in Supplementary Figs. and to be compared with the DD-vdW and DD results. It is seen in Supplementary 
Fig. 1^ that the mean-field simulation fails to reproduce the observed phase-shift again for both the DD-vdW and DD 
interactions. It is seen in Supplementary Fig. |12[ on the other hand, that the theory model simulations beyond mean- 
field with the DD-vdW interaction show good agreements simultaneously with the Ramsey contrasts and phase-shifts, 
whereas similar agreements have not been found with the DD interaction. Supplementary Figure shows the two- 
atom potentials that we have obtained by diagonalizing the ®^Rb Hamiltonian for two atoms with the DD interaction. 
It is seen from this figure that there is a mixing of states correlating to the 42 D 5 / 2 + 42 D 5/2 and 44 P 3 / 2 + 4 OF 7/2 
asymptotes due to DD coupling at interatomic distances shorter than ~ 2 pm. This distance is longer than the 
shortest interatomic distance accessible by our ps pulse excitation, indicating the validity of the DD-vdW interaction. 
The DD-vdW interaction, however, needs two fitting parameters C 3 and rc, whereas the pure vdW interaction needs 
only one fitting parameter Cq. We have, therefore, employed the pure vdW interaction in the main text to suppress 
the ambiguity of the fitting. 

Supplementary Figure shows the simulated Ramsey contrast and phase-shift at r = 500 ps as functions of the 
cutoff radius rg (the lower abscissa), at which we truncate the integration in Eq. (6) in the main text, and of an 
average number of interacting atoms within the volume V = ^(rg — r^) (the higher abscissa). They are simulated 
by the theory model with the continuum approximation for the vdW, DD, and DD-vdW interactions without 
anisotropies. The population of the 42 D 5/2 is set to ~3.3% in these simulations. The converged values agree well 
with the measured ones indicated by the dark-grey solid lines, each of which is the average over eight points around 
T = 500 ps in Figs. 6b or d in the main text. The light-grey shaded area represents one standard deviation of the 
average over those eight measured values. On the other hand, the Ramsey contrast and phase-shift simulated with the 
DD interaction neither converge nor agree with the experimental observations. The black solid line in Supplementary 
Fig.[5ji shows the Ramsey contrast | 3 ('r)| given by Eq. (7) in the main text with 7 (t) = 0, which gives the situation 
in which the Ramsey oscillations modulated by the Rydberg-Rydberg interactions are dephased completely. This 
black solid line thus represents the upper limit of the contrast decay and accordingly the lower limit of the number 
of atoms ~ 32 to reproduce the Ramsey contrast ~ 0.45 measured at r = 500 ps, irrespective of the potential curves. 
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